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Abstract 

We consider a cosmology in which the final stage of the Universe is neither accel- 
erating nor decelerating, but approaches an asymptotic state where the scale factor 
becomes a constant value. In order to achieve this, we first bring in a scale factor 
with the desired property and then determine the details of the energy contents as 
a result of the cosmological evolution equations. We show that such a scenario can 
be realized if we introduce a generalized quintom model which consists of a scalar 
field and a phantom with a negative cosmological constant term. The standard cold 
dark matter with Wm = is also introduced. This is possible basically due to the 
balance between the matter and the negative cosmological constant which tend to 
attract and scalar field and phantom which repel in the asymptotic region. The 
stability analysis shows that this asymptotic solution is classically stable. 
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1 Introduction 



Recent observational data [H [2] suggests that dark energy is believed to be responsible for 
the observed acceleration of the Universe [5]. Many theoretical proposals on the nature 
of dark energy have been put forward. Cosmological constant [H El [6l [TJ |8] is widely 
accepted as the most promising candidate for dark energy, but there is the problem of 
its extreme fine-tuning. There are many alternative approaches which do not contradict 
with observational data [S]. One possible approach is based on the fact that the matter 
with the equation of state parameter w < — 1 is not ruled out by observational data and 
moreover necessary to describe the current acceleration of the universe. Phantom [10] 
is an alternative form of dark energy with this property and various aspects have been 
discussed extensively [HI [T21 [131 E] • 

In spite of many theoretical attempts, at present the nature of the dark energy is still 
enigmatic. An alternative approach to dark energy could be to try the "reverse" way; 
Instead of guessing the energy density and then trying to determine the evolution of the 
universe, one first assumes some specific scale actor a{t), and then fix the energy density 
from the the equations of motion [151 IlHl CZl [E] • This type of approach was proposed 
before in which the focus has been primarily on the description of the present universe 
starting from the past. In this work, we will be interested in the future of the universe. 
So far, many results of the existing literatures on the future of the Universe is based on 
the consequences of acceleration [IHl EQl El] . 

However, the ignorance of the exact content of the dark energy leaves other choices 
of the scale factor as viable. In this paper, we will introduce a particular form of the 
scale factor which is consistent with the present accelerating universe. It represents the 
universe which does not expand forever, but approaches asymptotically a constant valuj^]. 
The motivation for considering such Universe is that there is a priori no reason why 
the acceleration of the universe should be the final state. The A— CDM suggests the 
acceleration as a final stage, but other contents of the dark energy could give different 
fate of the Universe. For example, in the phantom model the energy density may increase 
with time and the Universe ends up in big-rip singularity [HI [121 ES] in which the Hubble 
parameter and the curvature diverges in a finite time. Another speculation is that dark 
energy might dissipate with time and Universe contracts in a big-crunch [271 128] as in 

"'^The difference with the Einstein static universe is that in this case, the universe stays static or 
starts from a static state [22], whereas in our case the static state is approached asymptoticahy in the 
future. It can be shown that the Einstein static universe is stable against inhomogencous perturbations 
[251 [211 [53], even though it is unstable when the homogeneous perturbations arc performed. We only 
consider homogeneous perturbations in this work. 
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the cyclic model [291 EOl EI]. Here, we investigate Universe neither accelerating nor 
decelerating, but becoming asymptotically static with its size approaching a constant 
value, and ask what type of energy contents will yield such Universe. 

In order to implement the idea, we consider a generalized quintom model [5^ [5^ 
which consists of a scalar field and phantom, and also the cold dark matter which satisfies 
the continuity equation. The potentials of the scalar fields and other quantities defining 
the model are not fixed at this stage. With this, we bring in a scale factor which has 
the property that it is accelerating at present, but begins to decelerate after a fixed time, 
and eventually settles down to an asymptotic value. Then, the details of the theory are 
determined through the evolution equations. We show that such a scenario is possible due 
to the balance between the attractive matter and the negative cosmological constant, and 
the repulsive scalar field and phantom in the asymptotic region. The stability analysis 
shows that this asymptotic solution is classically stable. 

2 Phantom model 

Let us first consider phantom cosmology where gravity is coupled to a scalar field with 
a scalar coupling function and a potential. This model was considered before [16] where 
the possibility of phantom-non-phantom transition occurs in early time as well as at late 
time. In this paper, we go one step further and discuss the future of the Universe with 
such a scalar coupling function. 
The action is given by 

S = J d'xV^{^R-^u{^)V,ipV''^-Vi^)] + Sm, (2.1) 

where = SttGn is the gravitational coupling, u}{ip) can change continuously, and V{ip) 
is the potential. 5*^ is the matter action which we assume to be a cold dark matter with 
Wm = 0, and satisfies the continuity equation. The scalar field (p is a function of time 
alone. The evolution equations are given by 

— U = Pm, 

2 . 

— -i^H = + + pm, (2.2) 
where H = a/a is the Hubble parameter, the energy density p^ and the pressure are 

= ^£l'(<^)0'-\^, (2.3) 
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and the scalar field satisfies 

Co{^)Cp + ]^Co'{^)^^ + ?,HCo{^)^ + = . (2.4) 

There could be many expressions of scale factor which has the desired property of 
transition from acceleration to deceleration. In this work, we choose a special ansatz of 
the scale factor as follows; 

a(t) = aoo(l-Ae-^(*/*°)"), (2.5) 

where A, B are constants, to is the current age of the Universe, a^o is the asymptotic value, 
and n is arbitrary. This ansatz shows the behavior that the Universe accelerates until 
some time (given below) and then deceleration takes place until it reaches an asymptotic 
value. Therefore, it does not fit the whole history of the Universe, and we assume it is 
an effective scale factor which describes the Universe only after the late-time acceleration 
has begun, ttr — O.Sto- The transition from acceleration to deceleration will happen at 
the time satisfying the following conditions: 



The condition that the deceleration must take place at later time restricts n > 1/(1 — 5). 
With this scale factor, the Hubble parameter is given by 

_ 1 AEn(t/to)""^e"^(*/^°)" 
l-Ae-^(*Ao)" 

and the matter density is given as 



Pmit) 

With the ansatz as in [16] 



Po(l - Ae 



-B\3 



(1 - Ae^^(*/*o)")3 ■ 



^ = t, (2.7) 



the Einstein equations (12. 2 p can be rewritten as 

1 /„„o... .V.A 1 
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In order to determine the constant A and B, we demand 



a(to) 



aoo(l - Ae ) = ao 



1 ABne-^ 1 



to 1 - Ae-^ to ' 



Then, we get the following two relations 
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B 



n{m — 1) 



where m = — 



The free parameters are given by m and n; m determined the final size of the universe, 
and n is related with the transition time from acceleration to deceleration. They can be 
restricted by the observational data. The effective equation of state parameter is given 



and at the current time we confront the equation of state parameter w^p with experimental 
uncertainty (—1.14 < < —0.88) [2] which constrains the parameter n and m. The 
numerical analysis is given in Fig. 1 which shows that there exists allowed region within 
the experimental uncertainty. Among the possible values, we choose a point and present 
numerical result in Fig. 2. The result shows that at current time, the potential is dominant 
so that the Universe is accelerating. Also the scalar function is slightly negative and it 
is in the phantom phase. As it evolves, the deceleration begins to take over as the scalar 
function a)((/?) (kinetic energy) reaches its near peak value. Then, both the potential and 
the scalar function decrease rapidly and becomes negative. However, the matter density 
stays always positive, decreasing as 1/a^, and asymptotically it cancels exactly with the 
scalar contributions of the potential and the kinetic energy. This exact cancelation renders 
the Universe asymptotically static. 

In Fig. O we plot various quantities for values of (m, n) = (3.5, 2). 



Let us check the stability of the ansatz ( 12. 5 p and (12. 7p . In fact, this model( with one 
scalar field) could be unstable due to several reasons. The first one is that the effective 
equation of state parameter diverges, which occur in case of vanishing of denominator in 



by 



P. _ ^ - vjt) 

P<P ^ + V{t) ' 




3 STABILITY OF THE PHANTOM FIELD 
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Figure 1: The experimental data constraints to the {n,m) parameter space. The blue shaded 
region are allowed by experimental uncertainty —1.14 < w^p < —0.88. (We set the scale factor 
uq = 1, current time to = 1) ^^nd k = 1.) 

(12.81) . The other one is the instabiUty of the solution (12. 5p and (12. 7p due to singularity 
which occurs when crossing the point of the coefficient of kinetic term oj{(f) becoming 
zero as discussed in [16]. In order to avoid these problems, we may consider two scalar 
fields model like 

S = J rf4xv/=^|^(i?-2A)-^a;(0)V^0V'^0-^r/(x)V,xV^— 1^(0, + (3.9) 

where V{(j), x) ^-^e potential energy for which we assume that there is no direct coupling 
between field and x field, V{(f>,x) = ^<p{(p) + ^xM- A is a cosmological constant. For 
our purposes, we take the value of A to be negative. 

By varying the action, one can obtain the Einstein equations and the continuity equa- 
tion for CDM as follows; 

3 rr2 ^ 

—H = + p.^ + p,„ + — , 
2 • 

" P't' + Px + P'P + Px + P-^ ^ (3-10) 

= p^ + 3Hp^, (3.11) 

where energy density p^, py^ and pressure p^, are given by 

P4> = \^i.<P)<P'^ + y<l>^ Px = \vi.x)x^ + V^, 
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Figure 2: The left graphs are time evolution of the a(t) (solid curve) and 
H{t) (dash — dotted curve). The right graphs are time evolution of the V{t) (green solid curve), 
Pm{t) (dash — dotted curve), and w (red solid curve). (The parameter values are used for n = 2 
and m = 3.5. We set the scale factor ao = 1, current time to = and k = 1.) In these figures, 
the Universe accelerates until the time t* = 1.58. 
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(3.12) 



and equations of scalar fields and x as 

= uj{(j))4) + ^uj'{4>)<i)^ + 3Huj 

1 . 





d4> 



V{X)X + lv'{x)x^ + 5Hr]{x)x + 

2 dx 



(3.13) 



In the two scalar fields model, we can choose the corresponding coefficients u and rj 
to be non-vanishing everywhere. Then the instability due to divergence would not occur. 
For example, we may choose that u should be always positive and r] be always negative 



w(0) = + w(0)2, 

Vix) = ^(x) - A/a2 +cl;(x)2 



(3.14) 



where a is a non-zero constant. Since u and t] are positively- and negatively-defined, re- 
spectively, u gives ordinary scalar phase and rj phantom phase. The numerical evaluation 
is shown in Fig. |3l The corresponding potential energy V^(0) and V^{x) are obtained by 
solving the following differential equations that are just the scalar ffeld equations ( I3.13p 
with solution (p = t and x = t 



1 

— u 
2 



3Huj, 
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1 

"2^ 



(3.15) 



where the solutions and are generally dependent upon integration constant which 
can be defined in accordance with the potential energy of the one scalar fields model. 
The corresponding effective equation of state parameter and are given by 



Px 
Px 



(3.16) 



As was mentioned before, the effective equation of state parameter of the one scalar 
field model becomes singular at some point. In order to obtain the well defined effective 
equation of state parameter of the two scalar field model, we can choose the integration 
constant such that denominator in (13.161) should be always non-vanishing. This can be 
done and some numerical result is shown in Fig. |4] 
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Figure 3: The graphs are evolution of the coefficient of the kinetic term uj{(j)) and ??(x)- The 
dashed curve is coefficient of the kinetic term a) (93) of the one scalar field model. (We set the 
parameter value a = 1). 

To gain some insight into the property of the solution, we introducing the following 
variables |16] : 

X^ = 0,X^ = X, (3.17) 

and we consider the perturbations of variables , , H , and pm near the background 
solution. One can obtain the equations for the linear perturbations up to the first order 
as follows: 

^5X^ = (-— - 3h] 6X^ - 36H , 
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Figure 4: The left graphs are evolution of the potential energy of (red curve) and 
(blue curve). The dashed curve is evolution of the potential energy of the one scalar field 
model. The right graphs are evolution of the effective equation of state of field (/> (red curve) and 
field X (blue curve). (We choosing the initial condition V^(X) = 6 and V^x(l) ~ ^(1)— Kii(l)— A/k^ 
with A = —6 and = 1.) 



dt 



---3H] 6X^-36H. 
V 

io6X, 



3p,n6H - 3H6pm ■ 



(3.18) 



Before we carry out the numerical analysis, let us analyze the late time property analyt- 
ically. When time t becomes very large, the contribution of H ,u' , and t]' will become 
negligible. Therefore Eq.( ]3.18p reduces to 



d 
di 



( ^^'f> \ 
6X^ 



SH 

V ) 

where bj^ , i]^ and 
by 



M 



( ^^<t> \ 
bX^ 

5H 



M 



( 



V 






-^700 





-3 \ 
-3 


' Pm,OQ 



(3.19) 



-p„,oo / 
are values at late time. The eigenvalues of matrix M are given 



0, 0,± 




2 Pm,OQ 



(3.20) 



Since + ^700 + |Pm,oo — is negative (for the positive value m), there are no 

positive eigenvalues. Thus, the solution is stable in the late time. We mention that this 
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result can also be obtained by solving the differential equations Eq.f l3.19|) . and in the late 
time region, the solutions are 



SX^{t) = A^ + B cos(At) + C sin(At) , 

5X^{t) = A^ + B cos(At) + C sin(At) , 

SH{t) = ^(^-Ccos(At) +5sin(At)) , 

6pUt) = Ap^ + Poo {B cos(At) + C sin(At)) , (3.21) 



where A, A^, Ay, Ap^^^, B, and C are given by 



A^ 
B 




C = — -i, (3.22) 



with Ji, I2, I3, and I4 being integration constants. We find that the solutions has an 
oscillating behavior with finite amplitude at late times. A direct numerical analysis of the 
Eq. (I3.18P also confirms the stability. The results are shown in Fig. |5l It shows that the 
perturbation settles down to the one given in Eq. (13.211) very rapidly. The subsequent 
variation oscillates around the solution. 



4 Conclusion 

In this paper, we considered a generalized quintom model which consists of a scalar field 
and phantom, and also the cold dark matter which satisfies the continuity equation. The 
potential of the scalar fields and other quantities defining the model are determined by 
the condition that they are consistent with the scale factor which has the property that 
it is accelerating at present, but begins to decelerate after a fixed time, and eventually 
settles down to an asymptotic value. The details of the theory are determined through 
the evolution equations. We showed that the asymptotic scenario is possible due to the 
balance between the matter and the negative cosmological constant which tend to attract 
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Figure 5: The left graphs are time evolution of the 6X(f, (solid curve) and 5Xy- (dashed curve). 
The right graphs are time evolution of the 5H (solid curve) and 5pm (dashed curve). (The 
parameter values used are n = 2 , m = 3.5 and a = 1 and computation done by choosing the 
initial conditions 5X^ = Q.mi,5X^ = Qmi,5H = 0.000165, and dpm = 0.001.) 

and scalar field and phantom which repel in the asymptotic region. The stability analysis 
shows that this asymptotic solution is classically stable. 

The scale factor showing the behavior of approaching an asymptotic state is by no 
means unique. We could try other ansatz which exhibits similar property with the one 
investigated in this paper. Obviously, for different choices of the scale factor, the energy 
contents will be given differently. Nevertheless, the over-all feature of necessity of the 
phantom will remain, because in order to achieve the asymptotic final stage, the attractive 
force of the matter must balance the repulsive phantom. It remains to be seen whether 
the exact nature of the dark energy could support such an asymptotically static Universe. 
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